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ABSTRACT 

Among the transient black hole binary systems, 4U 1630-47 is one of the most 
active sources exhibiting outbursts every few hundred days, with every outburst last¬ 
ing typically around hundred days. During the 2002-2004 outburst the appearance of 
quasi-periodic oscillations (QPO) coincide with the onset of anomalous state. There 
are two distinct QPO states: namely the single QPO state with one QPO and the 
twin QPO state with two QPOs not related harmonically. The spectral features of 
this state are corroborated by a previous outburst in 1998. The evolution of the inner 
disc temperature and the inner disc radius suggest the possible onset of geometrically 
thicker, slim disc as a possible explanation of the energy spectral features. The other 
possibilities involving the Comptonizing cloud also exist that may explain the anoma¬ 
lous state. The two different QPO states exhibit different spectral features, and we 
provide a possible empirical physical scenario from the observation of the evolution of 
the spectral features. 

Key words: accretion, accretion discs - black hole physics - stars: individual: 4U 
1630-47 - X-rays: binaries - X-rays: individual: 4U 1630-47 


1 INTRODUCTION 


Since its first observations from Uhuru followed by Ariel V 
in the early seventie s, the X-ray binary source 4U 1630- 
47 (I Jones et al.l 1 19761 1 has exhibited reg ular outbursts wit h 
a gap of, typically, few hundred days (IPriedhorskvl [l98^ . 
The nature as well as the d urations of these outbursts can 
vary greatly (see for example lAbe et ^l2005l : lTomsick et al.l 
I 2 OO 5 I I. hence the study of spectral and the timing properties 
of the source can play a very important role in the fundamen¬ 
tal understanding of the accretion process for the transient 
X-ray binary systems in general. 

Neither the binary orbital period nor the na¬ 
ture or mass of the companion of 4U 1630-47 is 
know n (see lAugusteiin. Kuulkers. fc van KerkwiikI 
I 2 OOII. for a possible identification of the infrared coun¬ 
terpart). Nevertheless, based on the v arious obser- 


vations 

in the X-rays from 

Einstein dParmar et al.l 

T^. 

EXOSAT dParmar. 

Stella. & Whitd 

19861: 

Kuulkers 

van der Klis. & Parmar 

Il997bll. 

Ginqa 


ROSAT dParmar. Angelini. fc Whitel | l995lh Beppo -SAX 

llOosterbroek et al.lll998l l. ASCA i Parmar et al.lflQ^ 'l and 
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RXTE{ 

Kuulkers et al.l Il998l: iTomsick. Laoshov. & Kaared 

19981: 

Hiellming et al.l 1 19991: iTomsick & Kaared I 2 OOOI: 

Dieters et al. 2000l: 

Trudolvubov. Borozdin. & Priedhorskvl 

2 OOII: iKalemci et al. 

I 2 OO 4 II. iMcClintock & RemillardI d2003h 


classified the source as a ‘category A’ blac k hole candidate, 
suppo rting the earlier classification by iTanaka fc Lewii] 
(Il995h . 


Although the earlier observations of this source, includ¬ 
ing those till ROSAT studied mainly the outburst period¬ 
icity and the canonical state transitions, the power den¬ 
sity spectra (PDS) obtained from the lightcurves as ob- 
served by EX OS AT did suggest some ano malous behaviour 
dKuulkers. van der Klis. fc Parmaiill997bll. Thereafter, the 


spect ral observations from Beppo-SAX I Oosterbroek et aP 

Il998tl suggested the inadequacy of the optically thick 
and geometrically thin disc models under certain circum¬ 
stances. The RXTE/PCA observations of the outbursts in 
the e arly RXTE era, viz. 1998 outburst dHiellming et al.l 
I 1999 II brought the quasi-periodic oscillations (QPO) ob¬ 
served in th is source to the fore-front of study and 
classification (|D ieters_et_al, 2000||_ Tmsjclyfc KaarefII 2 OOC 1 I : 
iTrudolvubov. Borozdin, fc Priedhorskv 12Q01J). 


Among the nine outbursts exhibited by the source 4U 
1630-47 during the RXTE era (Fig.l[I]), the fifth one start¬ 
ing from -^MID 52525 to -^MID 53344 (2002 September 
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- 2004 December) 4U 1630-47 exhibited the lonRe st and 
brightest outburst phase observed till date. lAbe et alJ (l2005l i 
report the presence of three distinct spectral states includ¬ 
ing one anomalous state duri ng the outburs t s obse rved by 
RXTE from 1996 to 2004. iTomsick et al.l ll2005h report 
about the high amplitude variability, spectral states and 
transitions and also state that QP Os were observed dur- 
ing the 2002- 2004 outburst. Both lAbe et ^ (l2005l i and 
IXomsick et al.l ll2005ll suggest the presence of geometrically 
thick accre tion i.e. slim d i sc to be present during some ob¬ 
servations. iTomsick et al.l (l2005h also speculate that the un¬ 
physical low values of the inner accretion disc radius (as per 
the geometrically thin disc model) may be explained by the 
spectral hardening due to dominance of electron scattering 
in the inner region of the accretion disc giving rise to ob¬ 
served values of high temperature and luminosity. 

Despite the extensive coverage of this source in the 
literature, a comprehensive correlated study of the tim¬ 
ing features, namely the QPOs, and the spectral param¬ 
eters has not been carried out. In this paper we report 
the spectral properties of the source during the 2002-2004 
outburst viz-a-viz the QPO features observed. The re- 
sul t is corroborated by the measured parameter s reported 
by iTrudolvubov. Borozdin, fc Priedhorskvl (|200ll i from the 
1998 outburst. From the observational features of the spec¬ 
tral parameters and the QPO frequencies we conclude that 
the presence of low frequency QPOs occurs during a spectral 
state not classifiabl e by the present canonical black hole clas¬ 
sification scheme (jMcClintock fc R.emillardI l2003h and the 
possibility of the existence of a slim disc appears to be a 
viable alternative during the peak of at least a few of the 
outbursts of 4U 1630-47. Other alternatives may include 
non-standard Computerization in the innermost region of 
accretion flow. Nevertheless, this source is the perhaps the 
only one reported so far to demonstrate significantly a dif¬ 
ferent spectral state when the low-frequency QPOs appear. 


2 DATA AND ANALYSIS 


The All Sky Monitor fAS'M: [Levine et al.||l99 6|l aboard the 
RXT E satellite observatory ( Bradt. Rothschild, fc Swank! 
Il993h monitored the source continuously throughout its ex¬ 
istence, in the 1.5 - 12 keV range. In addition the nar¬ 
row field of view instrument, the Proportional Counter Ar¬ 
ray [RXTE/RCA) has covered the source extensively, espe¬ 
cially during the outburst phases. For the results reported 
in this paper, we have analysed 382 pointed observations 
of the RXTE/RCA (Fig. [2}, including 24 observations with 
the telescope pointing away from the source (offset ^ 1°), 
i.e. high-offset observations. The lightcurves were extracted 
from the event mode data for the 3-30 keV energy band and 
the complete band, at a time resolution of 2“® s. The energy 
spectra for all observations were extracted from standard 2 
mode data. To maintain uniformity only PCU 2 was used for 
the complete analysis, as it was the only common PCU avail¬ 
able for all observations. All the procedures of data filtering, 
background and deadtime corrections were strictly adhered 
to during the analysis process, including the cases of high- 
offset observations. The data reduction and analysis was car- 
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Figure 1. All the outbursts of the source 4U 1630-47 covered by 
the ASM during 1996-2011. The source shows a total of 9 out¬ 
bursts with a few outbursts exhibiting mult iple flares, the lo ngest 
being the outburst 5 during 2002-2004 ('see lAbe et al.|[2005h . 
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Figure 2. The Outburst 5 llAbe et*al.ll20n5h with multiple flaring 
episodes of 4U 1630-47 observed by the ASM. The rate is mea¬ 
sured in the units of Crab emission rate (75 counts s”^ of ASM). 
The arrows give the dates of PCA observations, the upward ar¬ 
rows represent observation for which no QPO is detected, and the 
downward arrows represent observations for which QPOs (one or 
more) are detected. 


ried out using HEAS OFldl, which c onsists of (chiefly) FTOOLS, 
XRONOS and XSPFC (lArnaudlll996l ). The power density spec¬ 
trum (PDS) were obtained from the event mode data for the 
complete energy band covered by the RXTE/RCA using the 
IDL based software package GHAT^. 

The PCA spectra were fitted for all the observations 
with a single three-component model consisting of mul¬ 
ticoloured disc blackbody dominating at lower energies 
{diskbb), a power law at higher energies and an isotropic re¬ 
flection component (reflect) of the power law emission from 


^ http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ 

^ package developed by T.M. Belloni at INAF - Osservatorio As- 
tronomico di Brera 
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the accretion disc (along with the interstellar absorption), 
as well as a two component model without the reflection 
component. Statistically good fits were obtained using the 
three component model for almost all the observations with 
the reduced x2 ^ 1-20 (with a systematic error ^ 1%) while 
fitting the model to PCA spectra in the 3-30 keV band us¬ 
ing the x^-minimization algorithm of the XSPEC package. 
For the complete set of observations considered, the best¬ 
fitting values of disc temperature Tin varies in the region 
0.7-3.2 keV, while the power law index F varies from 1.5 to 
a very steep value of 4.8. Although the possibility of any 
contamination from other sources in the energy range of our 
analysis is minimal des pite the location of the source be¬ 
ing in a crowded region dTomsick et al.ll2005li . we report the 
spectral parameters for the observations with 3-30 keV flux 
> 1.2 X 10~^erg cm~^s~^, just to ensure the removal of any 
effect on the spectral features due to contamination when 
the intrinsic flux from the source is weaker. 


The reflection component {reflect) was used to account 
for the broad excess around 8 — 10 keV for some observations 
in order to better co nstrain the inner disc te mperature Tn 
of the diskbb model JMakishima et aLlll98fill in XSPEC. Al¬ 
though this reflection component did not constrain the high 
values of Tin for certain observations, the x^ value was bet¬ 
ter constrained using this component for such observations, 
hence for the sake of consistency all the results reported in 
this paper are with the reflection component for all observa¬ 
tions. The lower limit of the column density parameter {Nr) 
was restricted to 6 x 10^^c m~^, which is the l owest mea¬ 
sured value for 4U 1630-47 l|Parmar et al.l IOOTII. Followin] 
the study of X-ray dips in the source ( Kuulkers et al.lll99i 
the value of 60° was used for the inclination angle, although, 
admittedly, there is huge uncertainty in this value as none 
of the binary parameters nor the companion star properties 
are ascertained. Similarly, the distance of the source was 
assumed to be lOfcpc. Consequently, the conclusions of 
this paper remain a qualitative statement based on the long 
term trends of the parameters as observed in the source. In 
case the actual distance or the angle of inclination is greater 
then, in either case, the value of Rin will be systematically 
greater. 


For certain high inner disc temperature Tn, the corre¬ 
sponding value of inner disc radius Rin is not always physical 
(Fig- El Section 5.1). For better constraint on the best-fitting 
values of these parameters, an attempt was made to use an 
empirical Comptonization model simpl which functions as a 
convolution with the diskbb component conver ting a fraction 
of th e input seed photons into a power law (ISteiner et al.l 
[ 200 ^. This model failed to fit the data statistically and 
the addit ion of additional Comptonizing component (namely 
comptt, lTitarchui3ll994li was needed to obtain an accept¬ 
able value. But even this model failed to constrain the 
values of Tn and Rin within physically acceptable limits. 
Also, a model consisting of only diskbb and comptt failed 
to constrain the values of Tn and Rin- Hence these models 
were discarded in favour of the model described in the pre¬ 
vious paragraphs to compare the current observations with 
the standard classification. 


In addition to using the GHATS, quite a few PDS (es¬ 
pecially for the large number of observations without any 


QPO ) were extrac ted using the ISI^ (iHouck fc Denicol^ 
l200d : jHouckI I2002H package with an add-on module called 
SITAto Of all the pointed observations analysed pertain¬ 
ing to the 2002-2004 outburst, statistically significant QPO 
was found for only 30 observations. All the PDS were fit¬ 
ted w ith single model consisting of peake d noise compo¬ 
nents JPsaltis. Belloni, fc van der Klislll999l l comprising of a 
zero-frequency centred Lorentzian (zfc) and four Lor entzians 
lH owa ^ [ 2000 ! : iBelloni. Psaltis. fc van der Klislf2002l l. as de¬ 
fined below: 

c Ro 

zfc = -- - - 

1 -I- vjvo 


and 


Lor = 


R?Qv<3 

vl T Q'^{v - uo)'^ 


where in the Lorentzian component, uq is the resonant fre¬ 
quency of the Lorentzian, Q{= vq/IS-Vq, where Aizo is the 
FWHM of the Lorentzian centred at vq) is the quality fac¬ 
tor and R is the normalization of the Lorentzian, such that 
the rms variability, rms =i?[0.5 — tan“^(—(5)/V]“^^^ (i.e. 
rms = R as Q —^ 0). The minimum criterion for the exis¬ 
tence a QPO is that the corresponding Lorentzian should 
have Q>1, while the preferable condition is that Q>2 and 
this latter condition was satisfied for most of the observa¬ 
tions. On many of the 30 observations with QPOs there exist 
multiple Lorentzians with Q ~ 2 or greater, but the centroid 
frequencies of the Lorentzians rule out the presence of har¬ 
monics or sub harmonics. The PDS was converted into a 
FITS file similar to energy spectra along with the creation 
of a dummy response matrix using the GHATS package and 
the x^ minimization technique of the XSPEC was utilized to 
obtain the QPO parameters. 


3 THE EXTENDED OUTBURST OF 4U 
1630-47 

During the RXTE era the source went into outburst phase 
multiple times, where some outbursts consisted of more 
than one flaring episode, including the 2002-2004 out¬ 
burst which i s referred to as outburst 5 (following the 
convention of lAbe et ^ l2005l f . A similar long outburst 
phase with more than one flaring episode probably oc¬ 
curred during 1988-1991 (total duration ~ 2.4 yr) as ob¬ 
served by Ginga, w hen the source flux r eached a peak 
value of ~ 0.6 Crab (iKuulkers et al.lll997ah . Furthermore, 
the source perhaps exhibited very bright outburst in 1977 
dChen. Shrader, fc Lividll997lf but the duration of the out¬ 
burst is not ascertained due to the sparse coverage of the 
source in that era of X-ray observatories. In comparison 
RXTE provides a continuous monitoring, via the ASM, as 
well as wide band spectral coverage coupled with high res¬ 
olution timing studies, via the PCA, from 1996 to 2011. 
Other than the 2002-2004 phase, the duration of the other 
outbursts is typically ~ 100 d (give and take few tens of 
days). Also, these typical short duration outbursts emit at a 


® http://space.mit.edu/CXC/ISIS 

http://space.mit.edu/CXC/analysis/SITAR 
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52500 52600 52700 52800 52900 53000 53100 53200 53300 53400 

Figure 3. The outburst no. 5 llAbe et al.l 1200,4) of 4U 1630- 
47 observed by the PCA. Bottom panel; 3-30 keV flux 
{xl0~^erg C7n~^s~^) as a function of days for which the PCA 
observations exist. Middle panel: the inner rad ius Rjn of the ac¬ 
cretio n disc obtained from the diskbb model llMakishima et ^ 
Il986tl . Top panel: the inner disc temperature Tin per the diskbb 
model. For all panels: the red marks indicate observations for 
which no QPO is present, the green and blue marks indicate ob¬ 
servations with single and twin QPOs, respectively. 

peak flux of ~ 0.4—0.5 Crab (|Abe et al.ll2005l : lTomsick et al.l 
|2005^ (see Fig. [TJ. 

The long outburst 5 comprised of five individual flaring 
episodes, each immediately following the previous one with a 
gap of few days (Fig.[2J, with the individual flaring episodes 
being referred to alphabetically A, B etc. The duration of 
flare 5A was from MJD 52525 to 52675, the flare 5B was dur¬ 
ing MJD 52686 to 52996, the flare 5C was during MJD 53004 
to 53061, the flare 5D was during MJD 53077 to 53251, the 
flare 5E was during MJD 53252 to 53344. The flares 5A and 
5B were evidently different from the other three flares as 
per the ASM data. The peak flux of these two flares reach a 
value ~ 0.6 and ~ 0.8 Crab, whereas the other three flares 
do not get brighter than ~ 0.3 — 0.5 Crab. Also, the peak 
phase of the flares 5A and 5B show large and fast fluctua¬ 
tions, whereas the other three show more stable behaviour. 
The arrows of Fig. [2] are the days for which the pointed ob¬ 
servation by PCA is available in the archives. The downward 
arrows indicate observations for which QPOs have been ob¬ 


served, while the upward arrows indicate observations with 
QPOs. 

In Fig. [3] bottom panel the 3-30 keV flux 
(x 10“®erg measured by the PCA during 

outburst 5 is plotted. The evolution of the inner disc radius 
and temperature, Rin and Tin, as obtained from the model 
component diskbb, are plotted in the middle and top panels, 
respectively. The data points of Fig. |3] are segregated as 
non-QPO (red mark) and QPO (green and blue marks) 
state. 

All the QPOs occur during the flares 5A and 5B, and 
that too only during the peak phase when the total flux as 
well as the spectral parameters show fast and drastic fluc¬ 
tuations. These fluctuations occur when, during the rising 
phase of the flare, the flux increases beyond ~ 0.3 —0.4 Crab 
flux. The fluctuations continue until there is a definite decay 
of the flare and the flux comes down below that same value, 
viz ~ 0.4 Crab flux. In between, during the peak phase fluc¬ 
tuations, the flux can reach a value less than ~ 0.4 Crab, 
and the parameter values of the inner disc is not discernible 
from the canonical black hole states. Hence, a more rigorous 
insight into the data is provided below to ascertain the be¬ 
haviour of the accretion disc of the low flux phases during 
the peak phase fluctuations. In this paper we have concerned 
our analyses only on averaged measurements of the different 
observational pointings. 


4 DIFFERENT CLASSES OF QPOS 

The observed PDS with QPOs are predominantly divided 
into two distinct classes, one with one QPO with the cen¬ 
troid frequency in the region of 11-14 Hz (henceforth re¬ 
ferred to as ‘single QPO state’), and the other with two non- 
harmonically related QPOs with the centroid of the lower 
frequency lying in the region of 4-6 Hz while the higher fre¬ 
quency can vary in the region of 7-18 Hz (henceforth referred 
to as ‘twin QPO state’). A sample each of the single QPO 
and the twin QPO is given in Fig. |4] In addition, another 
class of lowest frequency QPO with a periodicity around 1 
Hz is observed on a few occasions, and only once it appears 
alone in the absence of the other two classes of QPOs. Due 
to lack of statistically significant data, we will not include 
the lowest frequency QPO in our analysis. The details of 
the QPO parameters are given in table [T] The significance 
of the QPO was obtained from the errors calculated from 
the best-fitting statistics for the normalization parameter of 
the corresponding Lorentzian function. 

For the single QPO states, the peak rms as well as the 
significance is relatively greater for lower frequencies (the 
peak rms ^1% for frequencies < 12 Hz and the significance 
is highest for the lowest frequencies - table [TJ. For the twin 
QPO states both the peak rms as well as the significance 
shows large variation, but in general the peak rms of the 
lower frequency is equivalent or more than that of the higher 
frequency. The total flux measured from the source is more 
for twin QPO state except for the occasions when the higher 
frequency is in the range of 7 Hz (table [T|. 

In total there are 30 observations with QPO, 18 of 
them belong to the twin QPO state while 11 belong to 
the single QPO state, and once the lowest frequency QPO 
(« 0.84:Hz) appears alone. As per the prevalent ABC clas- 
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Table 1. The details of the observed QPO. The MJD (with the fractional part) denotes the start of the observation (there are occasions 
of multiple observations starting on the same day). The significance of the QPO is obtained from the best-fitting parameters of the PDS, 
see the text for details. 



MJD 

Obs. Id. 

QPO freq. 

FWHM 

rms-peak(%) 

sig. 

PCU2 count rate 

Single QPO 

52639.84286 

70417-01-09-01 

11.79±°'25 

2.39±i.09 

1 

2 

1254=1®! 


52639.92536 

70417-01-09-02 

11 81±0.43 

3 71±1.49 

1.1 

2 

1184±199 


52640.23910 

70417-01-09-03 

12.36±°-54 

4 41±1.64 

0.8 

4 

1278±46 


52793.38442 

80117-01-02-01 

12.53±oi2 

2 40±0.46 

1 

5 

172o±223 


52798.96622 

80117-01-04-00 

11.35±°03 

1 56±o.ii 

1.6 

16 

1757±194 


52800.84348 

80117-01-02-04 

12 4e±0.l3 

1 q3±0.42 

1 

5 

1804±213 


52810.09160 

80117-01-08-00 

11 4g±0.04 

0.94±o.ii 

1.9 

12 

1737±127 


52827.72771 

80117-01-11-01 

13.28±°'25 

1 14=1=0.61 

0.6 

2 

1Q29±149 


52834.50147 

80117-01-12-01 

11.60±oo6 

1 18±0.21 

1.7 

7 

1689±259 


52862.84906 

80117-01-16-00 

13.30±0'23 

1 95±o.81 

0.9 

3 

1791±153 


52868.08409 

80117-01-17-00 

13.10±°i® 

2,07±0.67 

0.9 

3 

1865±!™ 

Twin QPO 

52631.49685 

70417-01-08-00 

4.83±0-11 

0 91±0.32 

0.8 

4 

23g3±402 




13 27±0.34 

5 89±i.05 

0.7 

8 

2214±304 


52633.66694 

70417-01-08-01 


0.76±'^''*' 

0.6 

<2 




13.22±0'36 

5.60±2“ 

0.6 

2 

2164±217 


52636.79921 

70417-01-09-00 

4.62±“'^» 

2.14=^'^-'^^ 

0.5 

5 




12.54±0.16 

5 24±0.52 

0.7 

13 

2366*224 


52794.23388 

80117-01-02-02G 

4,71±'J'Jt> 

1.90±"-^' 

1.1 

12 




12.24*°°'^ 

3 44±o.20 

1 

23 

2410*!!® 


52794.91444 

80417-01-01-00 

4.77±'J.ui 

1.35±'^-^'‘ 

1.2 

7 




12.34±019 

4 96±o.57 

0.8 

12 

2424*!®® 


52795.28333 

80117-01-03-00G 

4,74±'J.'Jo 

2.08±"-^^ 

1.2 

9 




12.28±°°® 

3 26±0.24 

1 

18 

2494*!"“® 


52801.39948 

80117-01-05-00 

4.19±u-^u 

2 05±u.D» 

1.9 

4 




12.54±0.16 

3 70±o.46 

2.4 

11 

1383*!®® 


52802.79627 

80117-01-06-00 

5.61±“'^'=' 


2.9 

6 




7 36±o.07 

0 98±0.27 

2.4 

3 

1755*2'!® 


52802.86045 

80117-01-06-01 

5.67^"-'^“ 

5.00±'^-'*' 

2.1 

6 




7 82±0-04 

0.66±°13 

2.2 

6 

1690*303 


52804.23360 

80117-01-07-00 

5.72±“'"' 

4 72±u.a( 

2.1 

7 




7 56±o.05 

0 79±o.17 

1.8 

5 

2501*!®® 


52806.54943 

80117-01-07-01 

4 62±u.u4 

0 284=9.99 

1.1 

5 




15.15±°'3° 

6 29=11.35 

0.6 

4 

2581*!!2 


52808.24961 

80117-01-07-02 

4.53±'^''J^ 

0.594=912 

1 

7 




14 54±o.17 

6 81=10.53 

0.7 

18 

2538*!!® 


52812.26814 

80117-01-08-01 

4.31±“'^'= 

0.994=9 '>9 

0.7 

<2 




13.71±°'2"‘ 

5 11±0.79 

0.8 

7 

2275*!®® 


52819.22981 

80117-01-10-00 

4.53±'^''J'‘ 

0.304=9 12 

0.9 

3 




14.63±°-37 

6 46±i.56 

0.6 

3 

2366*®°® 


52820.82659 

80117-01-10-01 


0.324=9 13 

0.9 

3 




17.39±°'2® 

3.65±"‘-^2 

0.6 

3 

2521*!®° 


52824.29111 

80117-01-11-00G 


0 394=9.99 

1.1 

6 




14_41±0-22 

6 14±0.72 

0.7 

10 

2671*!®® 


52829.55105 

80117-01-11-02 


0 604=9.38 

0.7 

<2 




14.36±°'22 

4 18±i.03 

0.8 

4 

2476*!®! 


52864.86189 

80117-01-16-01 

4.08±"'“ 

3.734=1.12 

1 

2 




14 72±i.48 

12.28±o.77 

0.6 

3 


Misc. QPO 

52633.66694 

70417-01-08-01 

1 55±o.03 

0.18±o.08 

1.2 

3 

2214*®°® 


52802.79627 

80117-01-06-00 

1 39±o.09 

0.86=10.29 

3.1 

<2 

1383*1°® 


52802.86045 

80117-01-06-01 

1 39±o.09 

0 86=10.29 

1.8 

<2 

1755*2!® 


52804.23360 

80117-01-07-00 

1 27±0.05 

1 06=10.14 

2.4 

4 

1690*®°® 


52862.84906 

80117-01-16-00 

g 74±0.04 

0.23±o.ii 

2.1 

2 

1791*1®® 


52855.20911 

80117-01-15-01 

0 84±0.05 

0.35±o.15 

1.8 

3 

1729*20® 


Notes. The Miscellaneous QPOs have the lowest frequency and generally occur with other QPOs, except for MJD 52855.20911. The 

occurrences of single & twin QPOs are mutually exclusive. 
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Single Frequency QPO Twin Frequency QPO 



Frequency (Hz) 

Figure 4. Examples of the single and twin QPO states. 

sifica tion scheme for the QPOs (ICasella. Belloni. fc Stellal 
l2005h . the single QPO may be loosely classified as class C, 
despite the obvious deviation of the parameters from the 
class C values, simply because the centroid frequency rules 
out the A and B class. For the twin QPO state, the higher 
frequency deviates from C-type classification. Eventhough 
the lower frequency is nearer to the B-type, the QPO pa¬ 
rameters deviate significantly from the classification scheme. 


5 THE SPECTRAL BEHAVIOUR OF THE QPO 
STATES 

Most of the observations fall in the canonical black hole 
states, namely steep power law, thermal dominated and in- 
termediate while few fall in the hard state: as classified by 
iMcClintock fc RemillardI ll2003h . The anomalous states per¬ 
tain to the observations with QPO. Hence, for this source, 
the fundamental understanding of the evolution of the ac¬ 
cretion process needs to be constructed by segregating the 
spectral parameters of the QPO and non-QPO observations. 



Figure 5. The plot of inner disc radius versus inner disc 
temperature Tin for all the observations pertaining to flares 5A, 
5B, 5C, 5D, and 5E. The red marks represent the non-QPO states. 
The green marks represent the single QPO states and the blue 
marks represent the twin QPO states. The non-QPO observations 
obey the relationship oc T. (geometrically thin optically 
thick Shakura—Sunyaev disc). The QPO observations obey the 
relation Rin oc 


5.1 The Rin—Tin relation and its consequences 

The relation between the inner disc temperature. Tin and 
the inner disc radius Rin is plotted in Fig. [5] Pertain¬ 
ing to the five flares under study, flares 5A, 5B, 5C, 5D 
and 5E, the Rin-Tin plane has four branches. Two paral¬ 
lel branches (represented by red marks) correspond to the 
non-QPO observations from all the flares in the outburst. 
These two branches can be fitted by the empirical relation 
Rin oc . This is a clear si gnature of Shakura-Sunyaev 

disc llShakura fc SunvaevlIlQT^ ). The geometrically thin and 
optically thic k disc approximation is not violated for these 
observations (iKubota et al.l[l999l ). The other two parallel 
branches (represented by green and blue marks for QPO 
observations) correspond to high temperatures Tin and low 
inner-radius Rin and are best fit by a relation Rin « TiJ- 
This is the anomalous state which can not be modelled by 
a geometrically thin disc. Hence these observations provide 
the signature of the impending anomalous state which exists 
only along with the appearance of QPOs. 

Nevertheless, the values of the inner disc radius for 
the high temperature anomalous states, as obtained from 


the spectral fitting, are not physically acceptable. Although 
it is known that different physical mechanism causing the 
opacit y of the disc can give rise to a different spectral 
shape JShimura fc Takaharal ll9S^ ) which can cause an un- 
derestimation of the inner disc radius up to a factor of 5 
llMerloni. Fabian, fc Ros^ 120001 ) . the most likely cause of 
the breakdown of the diskbb model is t he formation of ge¬ 
omet rically thicker ‘slim’ accretion disc llAbramowicz et aP 
Il988l) . The Rin — Tn relation suggest that for a slim disc 
the emission at small radii can also produce a much flatter 
radial temperature p rofile due to more material in the disc 


radial temperature p r 

dWatarai et al.ll200g . i 


refer to fig. 4). 

Also, the departure from the standard multi-colour 
disc may be attributed to effects of electron scattering 
and corresponding changes in radi al temperature profile 
(iKubota fc Makishirnal 1 20041 ). Both iTomsick et al.l (l2005l ) 
and Abe et al.l ( 2005l Eeport the deviation o f the disc lumi¬ 
nosity from the relation Ldisc oc T^n- While iTomsick et al.l 
l|200flh speculates about electron scattering in the inner re¬ 
gion of accretion disc causing the high temperature and 
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hence deviation from the relation, lAbe et ^ (l2005li is 
more convinced that this deviation is due to the slim disc. 
The presence of QPO suggest a rotational symmetry in the 
process that gives rise to the quasi-periodic phenomena, 
whereas higher electron scattering would result in more wide 
band noise, with the exception of a formation of a shock 
dRao et al.ll2000l l. The so called ‘anomalous’ state may be 
classified as the QPO state, for this source, as it is with the 
appearance of the quasi-periodicity in the PDS the param¬ 
eters pertaining to the accretion disc show a discontinuous 
evolution. 

The two branches of the Ti„ — Rdn relation for the QPO 
states (Fig-O are resolved as two distinct QPO states, with 
the single QPO state giving rise to the branch with lower Tin 
and the twin QPO state giving rise to the parallel branch 
with higher Tn. The single QPO branch here includes the 
observation on MJD 52855. 


5.2 The powerlaw component 

The two distinct QPO states, single and twin, suggest that 
the dynamics of the accretion geometry is extremely var¬ 
ied. The discontinuity of the Tn — Rin behaviour indicate 
distinct accretion features hinting at different physical ori¬ 
gin of the single QPO and the twin QPO states, despite 
the similarity in the frequency and other parameters of the 
Lorentzian functions of these QPOs. The distinctness of the 
two QPO states is also seen in the plot of Tn as well as Rin 
w.r.t. the power law index P (Fig. [b]). A significant feature 
is that in the bottom panel the three green marks embed¬ 
ded in the non-QPO states correspond to the flare 5A. This 
shows that the for the single QPO state, the Tin forms a 
parallel track as a function of F in different flaring episodes. 
The corresponding points on the top panel {Rin) also shows 
a parallel track for the green marks corresponding to the 
same flare 5A. The twin QPO states do not form parallel 
tracks despite being present during both the flares. Thus, 
it is clear that the evolution of emission from the accretion 
disc is more complex and depends on more than one physi¬ 
cal origin, namely the accretion rate. The QPO states occur 
for a very narrow band of values of the powerlaw index, F. 



Figure 6. The plot of inner disc temperature Tin (top panel) 
and the inner disc radius Rin (bottom panel) versus the powerlaw 
index F for the QPO states. The red marks represent the non- 
QPO states. The green marks represent the single QPO states and 
the blue marks represent the twin QPO states. The three green 
marks forming a parallel branch correspond to flare 5A while the 
rest green marks correspond to flare 5B. The blue marks do not 
show any such demarcations. 


5.3 The fluxes of the spectral model components 

In Fig.[7]the scatter plot of fluxes of the model components, 
namely diskbb (bottom row) and power law (middle row) 
along with the ratio of the diskbb flux to the power law flux 
(top row) is plotted with the model parameters (Tin, Rin 
and F). The red plus marks represent the non-QPO state, 
the green circled dots represent the single QPO states and 
the blue circled dots represent the twin QPO states. It is 
evident that the diskbb flux for the different QPO states form 
parallel branches for the single and twin QPO states, and 
furthermore, the single QPO data points also form parallel 
branches for the flare 5A and flare 5B (the three green circled 
dots forming a parallel branch w.r.t. the rest correspond to 
flare 5A, while the rest of the green circled dots correspond 
to flare 5B). The twin QPO state, in contrast, doesn’t show 
any such demarcation despite being present during both the 
flares. This suggests that the emission from the thermal disc 
is very sensitive to the given local properties of flare in the 


case of the single QPO state in contrast to the twin QPO 
state. Furthermore, the power law flux shows a continuous 
evolution as a function of the inner accretion disc parameters 
as well as the power law index (middle row for Fig. [7)|. The 
power law index, F, doesn’t play any role in the emission 
from the disc blackbody (bottom right panel. Fig. [7]l, but 
even here the different flux values for the single QPO state 
of flares 5A and 5B form distinct parallel tracks. 

It is very significant from the Fig. that the power law 
flux is sensitive to the inner disk temperature, Tin, (Fig. m 
middle row) and not on the power law index, F. Thus, the 
power law flux is an excellent indicator of the presence of 
QPO state but not the power law index P. 

From Fig. [7] it is evident, from it’s relation with the 
diskbb flux as well as the power law flux, that the devia¬ 
tion of the Tin from its continuous evolution occurs during 
the QPO states. Furthermore, the QPO state is character¬ 
ized by less disc blackbody {diskbb) flux but more power 





















8 M Choudhury et al. 


10 


o cj 

2 o 

X ^ 

3 ^ 

EM 


0.1 


1.5x10- 


X 

E 


I 5x10-5 


10-8 

X 

E 

m 

m 5x10-5 


10-8 - 



+Non-QPO state 
oSingle-QPO state 
oTwin-QPO state 


O O ®00 

(330 


O ^ 

O O CfeO o oo" 

eP o 



Goo 


ogo 

qG 

®o 



^#+++++ ++t 

_^^^_I_^_I_^^^_I_L 

10 20 30 

R. (km) 


Figure 7. The plot of the fluxes of the different model components of the energy spectra versus the inner disc temperature Ti^ (left-hand 
column), the inner disc radius Ri„ (middle column) and the power law index P (right-hand column). The multicoloured disc blackbody 
(diskbb) flux is plotted on the bottom row, the power law flux on the middle row and the ratio of the two fluxes on the top row. The red 
marks represent the non-QPO states. The green circled dots represent the single QPO states, the blue circled dots represent the twin 
QPO states. The three green circles forming a parallel branch (distinctly for the diskbb flux and loosely for the total flux) correspond to 
flare 5A while the rest green circles correspond to flare 5B. 


law flux. Hence, the contribution of the standard geomet¬ 
rically thin disc appears to be minimal in the case of the 
QPOs seen in the source (top row, Fig. (Tj). One interesting 
inference from Fig. [7] is that the extreme steep power law 
observations do not correspond to the anomalous (QPO) 
state. Another important feature is that the disc blackbody 
flux evolves discontinuously as the source makes the transi¬ 
tion from the non-QPO state to single QPO state to twin 
QPO state and back; whereas, the power law flux under¬ 
goes continuous (and monotonic) change as the state transits 
across the three states defined by the presence and absence 
of QPOs. Thus, it may be conjectured that saturated Comp- 
tonization and the appearance of slim disc might happen to 
be two opposing features, each competing to establish itself 
at the other’s expense, and for this source the appearance 
of the QPO can be accepted as a signature of a geometrical 
evolution of the accretion disc. 


5.4 Correlation among the spectral parameters 
and fluxes 

To test for the correlation among the various model param¬ 
eters (namely Tin, Rin, F) and the fluxes of the model com¬ 
ponents {diskbb and the power law), we have used the Spear¬ 
man rank corr elation (SRC) c oefficient adapting the method 
prescribed bv iMacklirJ (Il982l) . For each correlation test be¬ 
tween two variables a set of three variables is chosen and 
the corresponding D-parameter is derived, which gives the 
confidence level, in terms of standard deviations, that the 
derived correlation between the first two variables is not due 
to the influence of the third parameter. If the D-parameter 
has a value > 111, then the possibility of the correlation of 
the two parameters can be ruled out as due to the third 
parameter. The important and significant cases from the re¬ 
sults of this correlation test is given in table [2] The tests 
were performed separately for the single and the twin QPO 
states (excluding MJD 52855). The obvious reason for the 
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separation of the two states is that the two states give rise 
to parallel branches in the parameter space, hence the same 
value of one parameter will give rise to discontinuous mul¬ 
tiple values of the other parameter, killing the opportunity 
of gleaning any useful information from the correlation test. 
Given the paucity of data points for the QPO states, the val¬ 
ues obtained for the SRC, and especially the D-parameter, is 
suggestive at its best, but nevertheless it provides important 
insight into physical scenario of the accretion process. 


5.^.1 Single QPO state 

From table [2] it is evident that all the reported correlations 
are very strong when only individual flaring episode is con¬ 
sidered (namely flare 5B; since flare 5A has only three data 
points its statistical parameters are not obtained, but a vi¬ 
sual inspection of Fig. [6] and [3 suggests a strong correlated 
behaviour of the various parameters) for the single QPO 
state. Since in flare 5B the Tin and Rin are perfectly anti 
correlated, the SRC and the D-parameter involving them 
and any other parameter have the same value with different 
sign. For the individual flare 5B, the only pair of parameters 
that give inferior correlation as compared to the combined 
flaring episodes is Tin and power law flux. This is a statisti¬ 
cal conundrum that requires more data, which is currently 
beyond the scope of the paper. 

While the anti correlation of Tn and Rin is obvious from 
the spectra (Fig. [3] and [SJ, the anti correlation of the Tn and 
the diskbb flux (or conversely, the correlation of the Rin and 
the diskbb flux) is a bit anti-intuitive, as it is expected that 
with the accretion disc extending further inwards the flux 
from the thermal disc component will increase in propor¬ 
tion, but just the opposite happens. In contrast, as the disc 
supposedly extends inwards, the power law flux increases 
(see values of Tn-power law flux and Rin-power law flux in 
table m also see Fig. [3). Furthermore, the power law flux 
has a very strong anti correlation with P, while the diskbb 
flux is very strongly correlated with F. This coincides with 
the fact that the power law index P is very strongly anti 
correlated with the inner disc temperature Tn and posi¬ 
tively correlated with the inner disc radius Rin, i.e. as the 
disc extends inwards the power law hardens and its emission 
increases. Conversely, the diskbb flux increases as the inner 
disc extends outwards and the power law index, P, steepens. 
Important point to be remembered is that the QPO states 
exist within a very narrow band of the values of F. Hence, in 
the QPO states, the power law index, P, is very sensitive to 
the accretion disc related properties and behaves in a very 
counter intuitive manner which is different from the normal 
behaviour seen in the non-QPO states (Fig. [6]and[3bottom 
right and middle right panels). 

The D-parameter suggests that the spectral parameters 
(viz. Tn, Rin, F) are all independent, while the correla¬ 
tion of the power law flux with the power law index, F, is 
most likely caused by the inner disc parameters. This sug¬ 
gest that the extent of the disc is more important a player 
compared to the presence of the Compton cloud, during the 
single QPO state. Additionally some significant change in 
the physical and/or geometrical state of the disc is being 
manifested in the inner disc parameters and their correla¬ 
tion with the model component fluxes. 


Table 2. The important and significant cases of the result of the 
Spearman partial correlation test among the model parameters 
and the fluxes of the model components, performed and presented 
separately for single and twin QPO states. 

Single QPO state (for flares 5A and 5B)(no. of data points: 11) 


Parameters 

(3’’'^ par) 

SRC coeff. 

Null prob. 

D-par 

^in ‘Rin 

(P) 

-0.88 

3 X 10“^ 

-4.83 

RinP 

(Tn) 

0.16 



Tn-T 

(Rin) 

0.23 



Tin". diskbb flux 

(Riu/T) 

-0.33 

~ 0.32 

>|1| 

Tin'.power law flux 

(Rin/T) 

0.95 

~1X10“® 

>|1| 

Rin'. diskbb flux 

(Tn/T) 

0.64 

~ 0.03 

>|1| 

Rin'.power law flux 

(Tn/T) 

-0.82 

~ 2 X 10“® 

>|1| 

T: diskbb flux 

(Tn/Rin) 

0.64 

~ 0.03 

>|1| 

r'.power law flux 

(Tn/Rin) 

-0.82 

~ 2 X 10“® 

>|1| 

Single QPO state for 

only flare 5E 

(no. of data 

points: 8) 


Parameters 

(3’’'* par) 

SRC coeff. 

Null prob. 

D-par 

Tin '‘Rin 

(P) 

-1.00 

0.00 

—oo 

RinP 

(Tn) 

0.93 

8 X 10“"* 


TnT 

(Rin) 

-0.93 

8 X 10-"“ 


Tin'. diskbb flux 

(P) 

-0.90 

2“® 

-1.08 

Tin'.power law flux 

(P) 

0.86 

6 X 10“® 

1.06 

Rin'. diskbb flux 

(P) 

0.90 

2“® 

1.08 

Rin'.power law flux 

(P) 

-0.86 

6 X 10“® 

-1.06 

F: diskbb flux 

(Tn/Rin) 

0.88 

3 X 10“® 

0.53 

F'.power law flux 

(Tn/Rin) 

-0.81 

0.01 

-0.14 


Since Tn and Rin are perfectly anti-correlated during the flare 


5B, the SRC and the D-parameter involving them have the same 
value with different sign (same as the SRC value). Also if either 
of Tn and Rin is one of the main parameters to be tested and 
the other is the third parameter, then D-parameter is undefined. 


Twin QPO state (no. 
Parameters 

of data points: 18) 

(3^^ par) SRC coeff. 

Null prob. 

D-par 

Tin •T.in 

(p) 

-0.91 

1 X lO-’’ 

-4.12 

RinP 

(Tn) 

0.83 

2 X 10“® 

2.41 

TnT 

(Rin) 

-0.74 

5 X 10-"“ 

0.35 

Tin'.diskbb flux 

(Rin/T) 

3 X 10-®5 


>|1| 

Tin'.power law flux 

(Rin/T) 

0.20 



Rin'.diskbb flux 

(Tn/T) 

0.80 

7 X 10-® 

>|1| 

Rin'.power law flux 

(Tn/T) 

0.22 



F: diskbb flux 

(Rin) 

0.62 

0.36 

0.30 

F: diskbb flux 

(Tn) 

0.62 

6 X 10-® 

1.03 

F'.power law flux 

(Rin/Tn) 

0.23 




5.4-8 Twin QPO state 

During the twin QPO state, the one parameter that is com¬ 
pletely different in values from the single QPO state is the 
inner disc temperature Tn (Fig. (3- In general, the data 
points in the scatter diagram of parameter values and the 
model component fluxes are more clustered than being corre¬ 
lated for the twin QPO state (table[3and Fig. [3, especially 
for the relation of flux values with Tn- Also, the clustering 
in different flaring episodes do not give rise to discernible 
branches in the scatter plots, nor are the correlation param¬ 
eters significantly altered if the analysis is done separately 
for the individual flaring episodes. But, since the cluster¬ 
ing of the data points happen in distinct regions, it is quite 
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possible that the QPO states have different physical origins. 
Yet, the general trend suggest somewhat similar evolution 
of the disc and the Compton cloud and their emission as to 
the single QPO state. 

The model parameters Tin, Rin and P show equally 
strong and independent correlation result as in the case of 
single QPO state. The significant difference of the SRC val¬ 
ues from that of the single QPO state is that the power law 
flux is not correlated to any of the model parameters, while 
the diskbb flux has weaker correlation compared to the sin¬ 
gle QPO state. The D-parameter suggest that the T-.diskBB 
flux correlation is due to the third parameter Rin (and in¬ 
dependent of Tin). 

Within the twin QPO state, the power law flux reaches 
a saturation w.r.t. Tin, while the diskbb flux falls with in¬ 
creasing Tin (Fig- O- But, in comparison with the single 
QPO state, the flux of both the components {diskbb and 
powerlaw) is higher for the twin QPO state. From the top 
row of Fig. [7] it is evident that during the twin QPO state 
the diskbb flux increases discontinuously. Hence, the twin 
QPO state, in a sense, depicts a more extreme version of 
the evolution of the accretion disc and the Compton cloud. 


6 FREQUENCY DEPENDENCE OF 
SPECTRAL PARAMETERS 

The QPO frequency dependence of the various parameters 
and flux values are essential for further insight to the ac¬ 
cretion system. In Fig. [S] the scatter diagram of the various 
spectral parameters {Tn, Rin, P, Powerlaw norm) are plot¬ 
ted with the centroid frequencies of the various QPOs. In 
Fig. [3 the scatter diagram of the flux values of the spectral 
model components are plotted with the frequencies. In con¬ 
cordance with the results mentioned in the previous section, 
the spectral parameters as well as the flux values cluster in 
different regions for the different frequencies of the QPOs. 
The SRC coefficient is computed for the correlation between 
the centroid frequency of the various QPOs and the spectral 
parameters as well as the flux values. The SRC is computed 
using any dummy 3'^'* parameter as in this case the computa¬ 
tion of the D-parameter is not helpful. The SRC coefficients 
are tabulated in table [S] 


6.1 Single QPO state 

The correlation of the frequency with the various parameters 
in this state is much stronger when the individual flare 5B is 
considered (table[3|. The frequency scales very strongly with 
the Rin and inversely with the Tn ■ Therefore, in the course 
of the evolution of the disc, the QPO frequency decreases 
as the disc extends inwards and subsequently the inner disc 
temperature Tn rises. The correlation of the frequency with 
power law index, F, is very strong, but weaker in compari¬ 
son to the correlation with the inner disk parameters. The 
frequency is also very strongly correlated to the diskbb flux 
and strongly anti correlated to the power law flux. 
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Figure 8. The various parameters of the model components that 
define the energy spectra for the QPO states versus the frequency 
of the QPO. Plus marks represent the single QPO state; filled 
triangles represent the higher frequency of twin QPO state; circles 
represent the lower frequency of the twin QPO state. 
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weak anti correlation (SRC=-0.56) between themselves. The 
behaviour of the higher frequency (7-18 Hz) is somewhat 
similar to the single QPO frequency, with the SRC coeffi¬ 
cient being less significant as the scatter of the parameters 
is more in the twin QPO state. Nonetheless, the frequency 
scales with the inner disc radius Rin and inversely so with 
the inner disc temperature Tn- But, the correlation of the 
frequency is highest with the power law index, P. This is in 
contrast to the single QPO state. This frequency is also cor¬ 
related with the diskbb flux, but with the power law flux the 
SRC coefficient value depicts a lack of correlation, hinting 
at a high scatter of this quantity w.r.t the frequency (but 
it is definitely not anti correlated, as was the case for the 
single QPO state). The D-parameter suggest this frequency 
is correlated to the Tn and the diskbb flux by the parameter 
Rin • 

The lower frequency (4-6 Hz) does not correlate with 
the model parameters (table[3]), although it forms a separate 
tight cluster in the scatter diagram (Fig.[8l. Intriguingly, the 
flux of the model components, however, are strongly anti¬ 
correlated to this frequency. 


6.2 Twin QPO state 

The twin QPO state show very distinct behavioural fea¬ 
tures of the two frequencies. The two frequencies show a 


7 EVOLUTION OF THE ACCRETION STATES 

Given the results obtained in the previous sections, it is in¬ 
teresting to note the evolution of the structure of the accre- 
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Table 3. The Spearman partial correlation (SRC) coefficient of the correlation between the QPO frequencies and the model parameters 
and the model component fluxes. The value of D-parameter given in brackets with the 3^*^ parameter when it is less than 1. 


Frequency type 

T- 

in 

Rin 

r 

diskbb flux 

power law flux 

Frequency: single QPO state 

-0.44 

0.68 

0.69 

0.93 

-0.34 

(both flares 5A and 5B) 

Frequency: single QPO state 

-0.95 

0.95 

0.83 

0.90 

-0.76 

(flare 5B only) 

Higher frequency: twin QPO state 

-0.80 {-0.55/Rin) 

0.85 

0.88 

0.70 (0.32/Rin) 

0.47 

Lower frequency: twin QPO state 

0.27 

-0.40 

-0.48 

-0.69 

-0.79 


Notes. The SRC coefficients are obtained using any dummy parameter when D-parameter is greater than 1. 
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Figure 9. The flux values of the model components that define 
the energy spectra for the QPO states versus the frequency of 
the QPO. Plus marks represent the single QPO state; filled tri¬ 
angles represent the higher frequency of twin QPO state; circles 
represent the lower frequency of the twin QPO state. 


tion disc. Starting from the quiescent state, the steady build 
up of both the flares 5A and 5B (Fig. result in the stan¬ 
dard Shakura-Sunyaev disc extending slowly inwards, along 
with the presence of the Comptonizing cloud (the power law 
component). The flares 5C, 5D, and 5E also behave in the 
same manner. But for flares 5A and 5B the rise of the flare 
does not stop with the luminosity peaking at ~ 0.4 Crab 
flux (like it does for flares 5C, 5D, and 5E), but continues to 
increase. This most likely leads to the QPO states when the 


applicability of the geometrically thin disc model becomes 
entirely infeasible. 

The temporal evolution may consist of the following 
stages. 1) The Shakura-Sunyaev disc extends inwards be¬ 
yond a certain threshold with the fall in and rise in Tin. 
The thin disc emission reduces and the power law emission 
(which can be called a non-thin-disc emission) rises. The 
onset of single QPO state happens and as the Rin falls the 
value of the single QPO frequency also falls while the Tn 
rises along with the hardening of the power law component, 
and the diskbb flux correspondingly falls while there is a rise 
in the power law flux. 2) The process may undergo a turn 
around and start evolving in the opposite direction, with a 
rise in Rin and a corresponding fall in Tin, in addition to 
all the other correlated changes in the parameters of the 
model and the QPO, resulting in the corresponding increase 
in the diskbb flux and decrease of the power law flux. 3) 
With further change in the accretion and/or disc parame¬ 
ters in the single QPO state, the onset of the twin QPO state 
occurs. The spectral parameters behave similarly, albeit at 
a definite higher value of Tn and lower value of Rin, sug¬ 
gesting that the accretion disc extends further inwards. The 
power law index P remains within the narrow limit analo¬ 
gous to the single QPO state. But the power law flux does 
not depend on the spectral parameters, while the diskbb flux 
now anti-correlates with P, the likely cause being the inner 
disc radius Rin- Here the higher frequency of the two QPOs 
behaves somewhat similarly to the frequency of the single 
QPO state, i.e. it scales with the Rin and inversely with 
Tin, but the correlation of the frequency is highest with P. 
It is the lower frequency that causes the system to behave in 
an inexplicable manner. The value of the frequency is anti¬ 
correlated with both the thin disc flux (diskbb) as well as 
the power law flux. 4) The system may fluctuate within the 
twin QPO state with corresponding changes in the parame¬ 
ters and fluxes, or it might transit down to the single QPO 
state. From the single QPO state it might again make a tran¬ 
sition to the twin QPO state or it may make the transition 
back to the canonical state as the flare dies down. 

Figs [2] and 1^ (bottom panel) show very fast fluctuation 
of the X-ray emission during the peak phase of the flares 
5A and 5B, and this feature is reflected in the spaced out 
data points of the QPO states in the Rin-Tn plane (Fig. 
EJ. This fluctuation of the flux (Figs [2] & [Sjl occur indepen¬ 
dent of any change in the mass accretion rate as the slope 
of the log-log diagram of Rin-Tn (Fig. El) remains constan t 
(IWatarai et aLlEoPOl : IWatarai, Mizuno, fc Mineshigel [20011 1. 
provided a transition from the geometrically thin to geo- 
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metrically slim disc has taken place. In this case the fluc¬ 
tuation of luminosity may suggest some instability in the 
accretion disc causing the drastic fluctuation of the Rin and 
Tin along with the onset of the QPOs, or else a fast physical 
mechanism is removing the mass from the inner disc close 
to the event horizon, resulting in the truncation of the disc 
and lowering of luminosity (Fig.Qbottom row). Eventhough 
the geometrically thin disc model is incapable of properly 
parametrizing the comparatively thicker slim disc parame¬ 
ters, but the total luminosity measured from the light curves 
(3-30 keV for PCA and 1.5-12 keV for ASM) is independent 
of the spectral model, suggesting that the fluctuation of the 
disc is a real phenomenon. 

This evolution may perhaps indicate that the geomet¬ 
rically thin disc component is actually incapable of giving 
proper explanation of the spectral evolution, and the power 
law component of the total flux may contain the emission 
actually originating from the thicker part of the slim disc. 
This may also explain the rise of the power law flux with 
the rise in Tin and fall in Rin as well as the diskbb flux. 
The onset of the geometrically slim disc may also explain 
the residue around 8-10 keV that is seen in the QPO states 
(see Section 2). 


7.1 Comparison with 1998 Outburst 


The 1998 outburst of 4U 1630-47 (outburst 2 as 
per lAbe et al.l l2005ll is well re ported in t h e lit- 
erature I Hjelinung et^ Il999l : fpieters et al.l l2000l : 
iTomsick fc Kaared l2nonl b The duration of this out¬ 
burst is around 100 d with one major flare with peak 
flux ~ 0.4 Crab lasting for about 60 d and a very 
small follow up flare with peak flux ~ 0.1 Crab flux, 
that lasted for about 40 d. Using the RXTE/PCA data 
iTrudolvubov. Borozdin, fc Priedhorskvl (l200ll ~) reported the 
presence of QPOs during the peak phase of this outburst 
which lasted for about 30 d. They also reported the 
best-fitting energy spectral parameters of the same (se e 
table 2 of ITrudolvubov. Borozdin, fc Priedhorskvl 1200 j l. 
The model used was essentially very similar to our model 
used here, the only difference being the addition of a 
cut-off power law to fit the 19-100 keV data obtained from 
i^ATE/HEXTE; the RXTE/PCA was used to obtain the 
spectra in the range of 2-20 keV. From their best-fitting 
parameters we have calculated the inner disc radius, Rin for 
the QPO-observations of the 1998 outburst, assuming the 
source to be at distance of 10 kpc and the binary system 
inclined at an angle of 60°. The resultant Tin — Rin plot 
for the 1998 outburst is plotted in Fig. 1101 where the data 
for the 2002-2004 outburst (erstwhile Fig. [5]) is also plotted 
for comparison. The three parallel lines representing the 
QPO states stand out in Fig. 1101 which suggest that the be¬ 
haviour has extremely close resemblance (albeit qualitative) 
with the behaviour expected in the presence of a slim dis c 
dWatarai et al.l[200ol : IWatarai. Mizuno. fc Mineshigel[200lll . 
During this outburst the presence of harmonically unrelated 
multiple Q POs in individual observations is reported (se e 
table 3 of ITrudolvubov. Borozdin, fc Priedhorskvl l200lf) . 
Hence, the outburst 2 might be very similar to the outburst 
5 (flares 5A and 5B) at a lower mass accretion rate. 

Interestingly, the follow up flares (5C, 5D, 5E, and also 
perhaps 2B) do not exhibit any quasi-periodic behaviour, 



Figure 10. The plot of inner disc radius Rin versus inner disc 
temperature Tin for all the QPO observations of the 1998 out¬ 
burst (peak phase)along with the 2002-2004 outburst (Fig. (SJ. 
The values of Rin and Tjn are obtained from the spectral p aram - 
eters reported by ITrudolvubov. Borozdin. &: Priedhorsk\l ll200lh 
(see the text for details). The red marks represent the QPO 
states of the 2002-2004 outburst; the green and blue marks rep¬ 
resent the QPO sates; the purple marks represent the QPO ob¬ 
servations of the 1998 outburst. The non-QPO observations obey 
the relationship Rin oc T. (geometrically thin optically thick 
Shakura—Sunyaev disc). The QPO observations obey the relation 
Rin oc T^ . 


while their energy spectral parameters also conform to the 
geometrically thin approximation. Perhaps outbursts with 
relatively strong mass accretion rate can only give rise to 
the quasi-periodic behaviour, which appears to be a strong 
indicator of the formation of slim disc. 


8 DISCUSSION 

The widely diverse nature of the emission from the Galactic 
black hole binaries require various classifications schemes to 
make sense of their radi ation mechanism and the underly- 
ing physical mechanism ijPone. Gierlihski. fc Kubotall200A 
and the references therein). The component based classifica¬ 
tion scheme classifies all the ‘possi ble’ spectral states of the 
black hole systems into five states llMcGlintock fc R.emillardI 
l2003lf . On the other hand the model independent clas- 
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(7-30 keV / 3-7 keV) Hardness Ratio 

Figure 11. The intensity-hardness ratio plot for the source 4U 
1630—47 during the outburst 5. The red plus marks represent the 
non-QPO state, the green crosses represent the single QPO state 
and the blue crosses represent the twin QPO state. 


sification scheme of I Homan fc Bellonil |200l) pays more 
stress on the transition between the states (Fig. 1111) 
and the related radio emission from the outflow from 
the system. The component based classification scheme 
llMcClintock fc RemillardI l2003h uses a multicoloured opti¬ 
cally thick geometrically thin accretion disc model (best 
parametrized by the simplistic diskbb model) with a com¬ 
bination of a power law and/or cutoff powerlaw. Hence the 
basic premise of the whole classihcation scheme is based on a 
geometrically thin disc. As a fallout, all the popular complex 
models used to fit the X-ray emission spectra are implicitly 
based on the assumptio n of the accretion disc being geomet¬ 
rically thin (see for eg. iK ubota et al.lll998l: lOier lihski et al.l 


I 1999 I : IZimmerman et a l.l l2005l:~ Gou e t al.| 201lh. Although 
the classification of Homan &: Bellonir i 2005ll is based on 
the colour-intensity evolution during the state transi¬ 
tion of the sources and hence independent of the spe¬ 
cific physical model used to describe the spectra, the 
physical interpretation of these transitions have an im¬ 
plied assumption of geometrically thi n disc llBellonil l2005l : 
iBelloni, Motta. fc Mufioz-Dariaj(201ll '). Comparatively, the 
attempt to ht observational data to the slim disc model has 
been sparse and few fnamelv IWatarai. Mizuno. fc Mineshigel 


The classifica tion of the low-frequency QPOs into 
three main types JWiinands. Homan, fc van der Klisl[l9^ : 

iHoman et al.ll200lli for the black hole binary system is quite 

established (ICasella. Belloni, fc Stellal 1200^ . This classifi¬ 
cation is helpful in providing greater depth to the under¬ 
standing of the energy spectral component based classifica¬ 
tions described above. This A-B-C classification scheme nor¬ 
mally requires the presence of one fundamental frequency 
along with its harm onic (or at times the sub harmonic) 
dCasella et akllioOdf l. The presence of multiple QPOs with 
frequencies not harmonically related in the source 4U 1630- 
47 adds to the complexity of classification of these states. 
But, the observation of simultaneous n on-harmonically re - 
lated QPOs is not unique to the source dMotta et al.ll20i^ . 
A complete understanding of all the QPO observations cov¬ 
ering the various types of outbursts of the source continu¬ 


ously monitored by the RXTE is needed to properly classify 
the QPOs obse r ved in this source. 

IPao et ^ d200nl '> do suggest that the low frequency 
QPO may occur due to the formation of shock in the tran¬ 
sition region of the geometrically thin but optically thick 
accretion flow and the low angular momentum optically 
thin accretion flow. Thus the possibility of high inner disc 
temperature occurring due to electron scattering can not 
be ruled out at the moment, with the seed photon of the 
Compton scattering originating at the optically thick geo¬ 
metrically thin disc and the Compton up scatter occurring 
in the optically thin cloud. Nevertheless, the failure of the 
Comptonzation models to ht the energy spectral data with 
physically meaningful parameters suggests that the Tin—Rin 
behavioural pattern in Fig. [S] and [TU] is qualitatively cor¬ 
rect. The similarity b etween the Tjn-Rin relationship and 
the model provided bv lWatarai et al] d200Cl[l provide the for¬ 
mation of slim disc as a very attractive physical option hith¬ 
erto not explored. Detailed physical mechanism of the thick 
accretion disc need to be modelled to explore the possible 
physical mechanisms that may g ive rise to the QPO phe¬ 
nomena (namely diskosiesmologv: IWagon^l2012l ~l. The pos- 
sibilit y of X-ray emissio ns from the jet during the anomalous 
state llTrigo et al.ll2013l l requires the explanation of the pos¬ 
sibility of the QPO phenomenon occurring in the outflow. 

The observed coincidence of the appearance of the 
QPOs and the onset of slim disc, with or without the outflow 
emitting in X-rays, opens up many new possibilities in the 
modelling of the accretion disc for this source in particular 
and other black hole binary systems in general. The need of 
the hour is to develop a generalized multicoloured accretion 
disc model that would enable a smooth transition from the 
geometrically thin to the geometrically ’thicker’, i.e. slim, 
disc, along with the emission from the out flowing jets. 


9 CONCLUSION 

In this paper we report the onset of an anomalous accretion 
state during at least two of the many outbursts exhibited by 
the black hole binary system 4U 1630-47. This anomalous 
state is very closely associated with the appearance of low 
frequency QPOs. The source of the QPO during the unusual 
state needs to be explored theoretically. One possibility is 
the formation of shock in the transition region of the ther- 
malized multicoloured geometrically thin and optically thick 
accretion disc and the Comptonizing plasma cloud which is 
geometrically thick and optically thin. The other option is 
the change in the physical and geometrical state of the ac¬ 
cretion disc and the thin disc transforming into a thick disc 
with increase in accretion rate. 

For this source the outburst may consist of more than 
one flaring events. During each flaring event, the mass accre¬ 
tion is perhaps the most legitimate cause of the peak value 
of the flux measured from the flare. From the two outbursts 
covered (1998 and 2002-2004), it appears that the onset of 
the slim disc, if it appears, happens if the flux approaches a 
value of ~ 0.4 Crab and beyond. Otherwise the geometry of 
the disc remains thin and the non-thermal emission is best 
explained to originate from the optically thin Comptonizing 
cloud, and/or from the shock region at the boundary of the 
accretion disc and the Compton cloud. In such a scenario the 
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formation of the slim disc and the existence of the Compton 
cloud might be competing events: but the last word in this 
can not be claimed from the given data reported here. The 
Tin-Rin behavioural patterns qualitatively suggest that for 
each flare the mass accretion perhaps does not change much 
but the slim disc phase sees a lot of fluctuation in terms of 
Rin as well as Tin- This change in the inner disc tempera¬ 
ture and the inner disc radius may be due to instability on 
the disc or the disappearance of matter from the inner part 
of the disc. 
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